QCD strings and the thermodynamics of the metastable phase of QCD at large N c 
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The thermodyanmics of a metastable hadronic phase of QCD at large Nc are related to properties 
of an effective QCD string. In particular, it is shown that in the large N c limit and near the 
maximum hadronic temperature, Th, the energy density and pressure of the metastable phase scale 
as£~ [T H - T)- (D ^' 6)/2 (for D± < 6) and P ~ (T H - T)- (D ^- 4)/2 (for D± < 4) where D± is 
the effective number of transverse dimensions of the string theory. It is shown, however, that for 
the thermodynamic quantities of interest the limits T — > Th and N c — > oo do not commute. The 
prospect of extracting D± via lattice simulations of the metastable hadronic phase at moderately 
large N c is discussed. 
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The modern incarnation of string theory as a funda- 
mental description of nature arose Phoenix-like out of 
the ashes of a failed attempt to describe strong interac- 
tion dynamics in terms of basic string degrees of freedom. 
With the advent of QCD, it rapidly became clear that a 
string model was not the correct way to describe fully the 
theory of strong interactions. However, string models did 
have many successes in describing important theoretical 
and phenomenological features of hadronic physics. On 
the phenomenological side, string models describe in a 
simple and natural way the linearly rising Regge trajecto- 
ries seen in the data[l| . They also give rise to a spectrum 
of hadrons which grow exponentially^ which is also seen 
in the datajJQ- On the theory side such models nat- 
urally have linear confinement. Intuitively this picture 
makes sense if color sources form flux tubes — as they ap- 
pear to on the lattice. Accordingly from the early days 
of QCD on, there has been a widespread belief that QCD 
in some sense becomes qualitatively equivalent to an ef- 
fective string theory, at least in a certain domain. Such 
an effective theory can be denoted as a "QCD string". 

This paper relates the nature of the QCD string the- 
ory to its thermodynamics at large N c . In the large N c 
limit QCD becomes a pure glue theory and has a first 
order phase transition between its hadronic and decon- 
fined phases 4] . The purpose of this paper is to show 
that at large N c the thermodynamics of the metastable 
hadron phase above this transition is determined by the 
nature of the string theory. In the first place regardless 
of the type of string theory there is a maximum tempera- 
ture Th above which the metastable hadronic phase can- 
not be maintained. Secondly, the power law behavior of 
the energy density and the free energy of the metastable 
hadronic phase is completely determined by the effective 
number of transverse dimensions of the QCD string. This 
is significant because lattice studies of the metastable 
phase of pure glue theory for relatively large N c (up to 
N c — 12) are now possible 0. 

The qualitative successes of the string picture can only 
be rigorously correct quantitatively in the large N c limit. 
For example, the linear Regge behavior and the expo- 
nentially growing density of hadrons both depend on 



the identification of hadron masses. However, the ex- 
cited hadrons are resonances and thus do not have well- 
defined masses. Fortunately, as N c — ► oo, the width of 
glueballs and mesons goes to zero (as A~ 2 and ATT 1 , 
respectively) 6], the ambiguity in the value of the mass 
vanishes, and the string theory prediction becomes well 
defined. Similarly, linear confinement only can occur in 
QCD if for some reason flux tubes are unable to break. In 
the large N c limit of QCD, a simple diagrammatic argu- 
ment shows that the formation of each quark-antiquark 
pair is suppressed by a power of 1/N C and thus at large 
N c , flux tubes do not break 6] . Therefore, the simple 
consequences ascribed to a QCD string theory can only 
accurately reflect the underlying theory in the large N c 
limit. 

Similarly, one expects a string-like description to be 
increasingly more accurate as the excitation energy in- 
creases. If the QCD string theory is simulating a flux 
tube, then one expects the dynamics of the tube moving 
collectively in the manner of a string to arise when the 
characteristic length of the flux tube is much longer than 
the typical width of the tube. Such excitations have large 
masses. 

For simplicity, the present analysis will be restricted 
to pure gauge theory. This restriction is motivated both 
by the focus on the large N c where quark loop effects 
are suppressed by factors of 1/N C) and by the desire to 
compare the analysis with lattice simulations which for 
relative large N c have only been done for pure glue||. 
The inclusion of quark loops leads to possible technical 
complications in the analysis and will be neglected here 
in the interest of clarity. 

Consider the behavior of large N c QCD at zero chemi- 
cal potential in the hadronic phase. As argued in ref. Q 
such a system is well described by a non-interacting gas 
of hadrons (which for our pure glue system means glue- 
balls). The hadrons are non-interacting since by standard 
A c counting agruments 6j glueball-glueball interactions 
scale with A c according to 

v ss ~ W? ■ (i) 

This scaling implies that if the number densities of the 
hadrons are of order N®, then the energy is dominated by 
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the mass of the hadrons also of order AT? and interactions 
may be neglected. 

A framework for describing such a non-interacting gas 
of hadrons has been developed for phenomenological rea- 
sons |8j and will be used here. Consider a gas of non- 
interacting glueballs where the k th species of glueball has 
mass m k and spin J k . It is straightforward to see that 
that the energy density e k and pressure associated with 
each spin state of the boson is given by 
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where K-x is a modified Bessel function and the ar- 
rows indicate asymptotic behavior of these functions as 
m/T — ► 00. The spectral density of hadrons is defined as 
p{m) = ~^2 k {2J k + 1) 8{m — m k ). The asymptotic contri- 
bution to energy density and the pressure in the hadronic 
phase are given by 



£{T) 



dm p(m)e(m,T) + 0(1/N 2 ) 



(3) 



where the (D(l/N 2 ) correction arises from the interaction 
energy between the glueballs. The specific heat can be 
obtained by differentiating the energy density with re- 
spect to the temperature. Generically both the energy 
density and the hadron number density are of order N®. 
At low temperatures and large N c , the energy density is 
dominated by the mass of the hadrons rather than their 
interaction energy; QCD can then be described increas- 
ingly accurately in terms of non-interacting hadrons as 
N c increases. 

However, at very high temperatures QCD is described 
increasingly accurately as a gas of weakly interacting 
quarks and gluons. At large N c this is dominated by 
the gluons and the energy density scales as N 2 . More- 
over, from the behavior of the (3 function it is apparent 
that the characteristic temperature at which the thermo- 



dynamics is well described in terms of weakly interact- 
ing quarks and gluons is of order N®. Thus, there are 
two temperature domains at order A,?: a low tempera- 
ture hadronic one with energy densities of order A° and 
a quark-gluon phase with energy densities of order N 2 . 
For infinite N c such a situation is only possible if there 
is a discontinuity: below some T c , £/N 2 = while above 
T c £/N 2 is nonzero. In fact, for pure glue theories we 
know from lattice studies that this transition is first or- 
der. Moreover, the latent heat of the transition appears 
to scale numerically with N 2 as expected from this simple 
analysis 0. 

Suppose that the high-lying spectrum of hadrons in 
large N c QCD is, in fact, describable as a hadronic string. 
The behavior of the density of hadronic states is known 
for asymptotically large values of M and is given by 



p(m) -> pasy(ra) = A (m/T H ) 



-2B e m/T H 



(4) 



where A, Th and B are parameters which are indepen- 
dent of N c . The value of the Tjj depends on the numerical 
value of the string tension a. It also depends on the type 
of string theory and, in particular, the world sheet cen- 
tral charge, c: T H = The keyto the present analysis 
is the parameter B is universal^ (i.e., independent of 
the type of string theory) and depends only on, D±; 



B={D± + 3)/4 . 



(5) 



The parameter T# is the well-known Hagedorn tem- 
perature associated with the spectrum 10]. In the 
hadronic gas model, the energy density in Eq. © diverges 
for any T > Th- With the advent of QCD it became 
apparent that the Hagedorn temperature was the upper 
bound for hadron matter [TH]; the system is necessarily in 
a deconfined phase above Th- As we now know, large N c 
QCD enters a deconfined phase before Th is reached^. 
However, Th presumably represents an upper bound for 
the temperature of a metastable hadronic phase Q. 

It was realized long ago that the power law behavior 
of thermodynamic quantities in the neighborhood of Th 
are determined by the power law of the prefactor in the 
Hagedorn expression; that is, by the coefficient i?[l^]. 
Evaluating the integrals for the asymptotic contributions 
to the energy density and pressure using the Hagedorn 
spectrum obtained from QCD strings yields: 



£ asy = (27r)- 3 / 2 AT( 9 -^/2 T 7 J 2 x 



(T H -T)(^- 6 )/ 2 r(^^) ifZ? ± <6 
- log {{T h -T)/Th)) if £>± = 6 
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P asy = (2n)- 3 / 2 AT( g - D ^/ 2 T 5 J 2 x 



Since the asymptotic contributions to £ and P diverge 
in the neighborhood of Th, they will dominate the full 
expressions which in turn will diverge in the same man- 
ner. The nonanalytic behavior of the energy density in 
the vicinity of Th immediately implies nonanalyticity for 
the specific heat. Quite generally one sees that 

a = 4-D ± /2, (7) 

where a is analogous to the usual critical exponent for 
the specific heat (c v ~ (T H - T) a ). 

Suppose, that the QCD string has D± < 6. The ther- 
modynamic behavior of such a system is qualitatively 
quite different from what one sees in typical statisti- 
cal mechanics models — at least at large iV c . For typical 
models the metastable phase can have a maximum tem- 
perature. Such a maximum temperature is a spinodal 
point — the point at which local fluctuations destabilize 
the phase. The energy density remains finite as the spin- 
odal point is approached — although the specific heat will 
diverge. In this sense, a spinodal point behaves essen- 
tially as a second phase order transition. If the QCD 
string has D± < 6, there is also a maximum tempera- 
ture for the metastable phase, namely, Tjj, as there is 
with familiar statistical mechanical models. However, in 
sharp contrast to the typical spinodal point, the energy 
density of the QCD string theory diverges as Tjj is ap- 
proached; for the case of D± = 2 one sees from Eq. (JHJ) 
that the power law behavior for the energy density of 
the metastable hadronic phase goes as £ ~ (T# — T)~ 2 . 
This inability to reach the maximum temperature in the 
hadronic phase due to its energetic inaccessibility is the 
modern incarnation of Hagedorn's old idea (although now 
relegated to a metastable phase and large N c ) and oc- 
curs for D± < 6. The critical exponent a encodes the 
difference between typical statistical systems for which 
the finiteness of the energy at the transition requires 
< a < 1 and the present case for which a is greater 
than unity (provided D± < 6). For example in the simple 
Nambu-Goto string (D± = 2) a = 3 . 

Lattice simulations showing clear evidence of the en- 
ergy density having such a power law divergent behavior 
for the energy density at large N c would provide strong 
numerical support for the notion of the QCD string. 
Similarly a numerical determination that a — 3 for the 
metastable phase would go a long way toward establish- 
ing that the effective QCD string has D± = 2 as in the 
Nambu-Goto string. 

Perhaps more exciting is the possibility that numeri- 
cal evidence could show that a / 3 near Tjj indicating 
the effective transverse dimension of the QCD string is 
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different from two. The effective transverse dimension 
of the QCD string could, in principle, be different from 
the number of transverse space-time dimensions since the 
effective number of transverse dimensions characterizes 
the number of degrees of freedom of the string. One can 
imagine scenarios in which the dynamics gives a larger 
number of degrees of freedom than in the simplest string 
picture. As a practical matter, properties of the QCD 
string can be extracted from the behavior of the system 
near the endpoint of the metastable phase only if D± < 6. 
For D± > 6, the energy density does not diverge at Tjj 
and it becomes difficult to disentangle the Hagedorn be- 
havior of string theory from the critical behavior seen in 
many statistical mechanical systems. 

An important caveat is in order. The preceding anal- 
ysis makes important use of the large N c limit. One 
might hope that the qualitative behavior of power law 
divergence encoded in Eq JSJ remains valid at large but 
finite N c (perhaps with the powers altered by 1/N C cor- 
rections). Unfortunately, for D±_ < 4 simple thermody- 
namics indicates that this cannot be the case. The pres- 
sure is the negative of the free energy density. By defi- 
nition in a metastable phase there exists another phase 
with lower free energy (and thus higher pressure). In this 
case it is the deconfined phase (with pressure ~ N 2 ). 
Thus, the expression given for the pressure in Eq. © 
breaks down when T gets sufficiently close to Th so that 
pressures are of order N 2 . This breakdown can be un- 
derstood physically — the analysis is based on negligible 
interactions between glueballs. However, as T — > Th 
the glueball density diverges. It is easy to show that the 
number density of glueballs scales in the same way as the 
pressure (at large N c ) and thus one expects that when 
the pressures reach order N 2 so do the number densities; 
Eq. then implies that the interaction energy does as 
well. Thus at that point the interaction energy is com- 
parable to the mass energy and analysis breaks down. 

It is clear that the temperature regime in which the 
expressions in Eq. JfjJ) dominate, the energy and the 
pressure is limited. For low temperatures, the asymp- 
totic form does not dominate, while very close to T c 
hadron-hadron interactions invalidate the analysis (pro- 
vided D± < 4). Evidently the domain of validity is 

f CaiVc" 3 ^ if,D ± <4 
ciT H «T ff -T«T ff x<^ (8) 

[ e~ C2N o if Dj_ = 4 

where ci and Ci are constants of order 7V^. Thus for the 
simple scalar string with D± = 2, the range of validity 
ciT H « T H - T « T H /N 2 . 
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Note that as N c — > oo the upper end of the region of 
validity approaches Th- However, for any finite value of 
N c the expression ultimately must break down as T ap- 
proaches Tjj- Thus the critical exponents for the infinite 
N c theory will differ substantially from the critical expo- 
nents for large but finite N c : the limit of the critical expo- 
nent as Nc — > oo will not equal the critical exponent for 
infinite iV c . Recall the critical exponent associated with 
some quantity f(T) that diverge at the critical point is 

defined as a limit: j f = lim T ^ Tc ((T c - T) dl ° s j^ {T)) y 

However, for the quantities of interest here the limits 
N c — ► oo and T — ► Th are not uniform. 

The nonuniformity of limits means that in extracting 
the power law behavior to learn about hadronic strings, 
one considers the N c — > oo limit before the T — > Th limit. 
Care must be taken in extracting the power law behavior 
of the energy density that one is working in the regime 
of validity of Eq. @, i.e., for temperatures not too close 
to Th- Since in practice it is very difficult numerically 
to study a metastable phase very near its endpoint, it 
may not be too difficult to do this provided N c is large 
enough. 

The study of iV c scaling behavior makes manifest the 
fact the range of validity of Eq. © is limited at the top 
end. This helps resolve a very old problem with Hage- 
dorn's original idea of a limiting temperature. The issue 
of whether hadron-hadron interactions could fundamen- 
tally alter the thermodynamics was recognized as impor- 
tant very early by Huang and Wcinberg|l2|. However, at 
the time there were no techniques to study the issue in 
a controlled manner. In practice it only could be han- 
dled in an ad hoc fashion by neglecting the interactions 
under the assumption that such interactions did not fun- 
damentally alter the physical density of states. The N c 
scaling arguments for the pressure show that this is not 
valid for the real world of N c — 3 or, indeed, for any 
finite N c — there will necessarily be a regime in which the 
hadron-hadron interactions must cause the growth in the 
pressure to saturate as T approaches Th while the ther- 
modynamics of a non-interacting gas of hadrons with a 
Hagedorn spectrum will have a diverging pressure. 

The preceding analysis provides no real insight into 
the behavior of the system at finite (but large) N c in 
the immediate vicinity of Th- A highly plausible sce- 
nario is that there still is a maximum hadronic temper- 
ature at a spinodal point T c which is very close to Th- 
((T h - T C )/T H ) ~ N^. As T c is approached glueball- 
glueball interactions become strong and this drives an 
instability leading to the existence of a spinodal point. 
As was argued earlier, for the case of D± < 4 the break- 
down in the hadron gas analysis occurs at temperatures 
where the energy of glueball-glueball interactions become 
comparable to the energies of the free particles; at this 
point, the dynamics changes qualitatively. Thus, when 
pressure of order N% are reached, nothing protects the 



system from developing an instability. 

The relationship of the nature of the QCD string ther- 
modynamics of the metastable phase of large N c pure 
Yang-Mills theory would be of comparatively little im- 
portance if there were no independent way to study this 
thermodynamics. However, lattice simulations of the 
metastable phase provide the possibility of direct calcu- 
lations of the thermodynamics from QCD. Simulations 
of the metastable phase at rather large N c have already 
been done[]|. These studies were aimed at identifying 
the maximum hadronic temperature (identified as the 
Hagedorn temperature) which was done by extrapolat- 
ing calculated values of the string tension down to the 
point at which it vanishes. If similar calculations were 
done for thermodynamic quantities such as the pressure, 
the energy density or the specific heat, it may be possible 
to get ab initio evidence of the nature of the QCD string. 

Note added after publication. After this paper was 
published, the author became aware of an important 
early paper by Charles Thorn on the subject of the QCD 
phase transition at large N c and its relation to the Hage- 
dorn spectrum 0. Reference concludes among other 
things that QCD must have a phase transition at infi- 
nite N c and that the Hagedorn temperature is bounded 
from below by the critical temperature for the transi- 
tion. Given the pioneering nature of this paper, and the 
fact that it is not as widely known as it should be, it is 
important to cite it here. 

The author is grateful to S. Nussinov for useful com- 
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